基于樟子松树轮晚材最大密度重建的大兴安岭北部7-8月平均气温数据集（1781-2013）

李明启^1*，兰 宇²

1. 中国科学院地理科学与资源研究所，中国科学院陆地表层格局与模拟院重点实验室，北京 100101；
2. 广东省气象台，广州 510640

摘 要：本文利用2013年9月份采自大兴安岭北部呼中国家自然保护区高海拔树线位置（51.79°N, 123.08°E, 海拔900 m）的樟子松（Pinus sylvestris var. mongolica Litv）树轮样芯（直径10 mm），使用瑞士生产的DENDRO 2003树轮分析系统获取了树轮密度数据，建立了树轮晚材最大密度（maximum latewood density，MXD）年表，并与漠河气象站气象数据进行了相关分析。结果表明MXD年表与该区7-8月平均气温相关最高。因此，利用MXD年表重建了1781-2013年大兴安岭北部地区7-8月的平均气温，重建方程的方差解释量为31.1%。数据集包括：（1）采样点的地理位置数据；（2）树轮MXD标准年表；（3）重建的1781-2013年7-8月的平均气温，及11年滑动平均值；（4）39根样芯的树轮晚材最大密度测量值统计特征。数据存储为.shp和.xlsx格式，由9个数据文件组成，数据量为27.2 KB（压缩为1个文件24 KB）。

关键词：树轮；晚材最大密度；大兴安岭北部；气温重建

DOI: https://doi.org/10.3974/geodp.2022.03.09

CSTR: https://cstr.escience.org.cn/CSTR:20146.14.2022.03.09

数据可用性声明：本文关联实体数据集已在《全球变化数据仓储电子杂志（中英文）》出版，可获取：
https://doi.org/10.3974/geodb.2022.04.02.V1 或 https://cstr.escience.org.cn/CSTR:20146.11.2022.04.02.V1.

1 前言

树木年轮具有定年准确、分辨率高、分布广泛、复本量大等特点，成为研究过去气候变化的重要代用指标[1]。在百年至千年尺度上区域气温[2]、降水/干湿变化[3-5]重建中起着重要作用。在用于气候重建的树轮数据中，树轮晚材最大密度（MXD）能很好地揭示生长季或者生长季末期温度变化[6-9]。东北地区是我国三大林区之一，森林面积达3,000多万hm²[10]，同时还是我国主要的粮食产区，2018年谷物产量占全国的20.26%[11]。温度是影响农作物和森林生长的主要因素，因此，研究我国东北地区气温变化历史，寻求其变化规律，对指导当地林业和农业生产具有重要意义。

在大兴安岭地区，树轮晚材最大密度已经被用来研究过去气温变化[12]，但是相对研究还较少，并且重建时段不到200年。本文拟利用采自大兴安岭北部地区的樟子松树轮晚材最大密度数据（图1和表1），重建该区过去233年气温变化历史，为预测该区未来气候变化提供基础资料，同时也为指导东北地区林业和农业生产提供科学依据。

图1 采样点位置图

下载

2 数据集元数据简介

《基于樟子松树轮晚材最大密度重建的大兴安岭北部7-8月平均气温数据集（1781-2013）》[13]的名称、作者、地理区域、数据年代、时间分辨率、数据集组成、数据出版与共享服务平台、数据共享政策等信息见表2。

3 数据研发方法

3.1 年表建立及气温重建方法

年表建立：用ARSTAN软件建立晚材最大密度年表，选择80年样条函数对每一条序列进行生长趋势拟合[15]，并用原始序列除以拟合趋势得到每条序列的标准化指数序列，最后将所有序列利用双权重平均法合成树轮晚材最大密度年表[16]。

相关分析：为了确定树木生长的限制因子，我们对树轮晚材最大密度年表与漠河气象站1959-2013年气候要素（月平均气温和月降水量）做了Pearson相关分析。结果表明：当年7-8月平均气温是树木晚材最大密度形成的限制因子。

相关系数计算公式：

(1)

式中，x和y分别为两个变量，δ为标准差。

气温重建：根据相关分析结果，利用一元线性回归方法建立大兴安岭北部地区樟子松树轮最大晚材密度年表（MXD）与7-8月平均气温（T_mean7-8）的转换方程，然后重建区域气温。

回归分析方程：

T_mean7?8 = aMXD + b (2)

式中，T_mean7-8是7-8月平均气温，MXD是树轮最大晚材密度，a和b为常量。

3.2 技术路线

我们于2013年9月份在大兴安岭北部东侧呼中国家自然保护区内采集了28棵树61根樟子松树轮样芯，样点位于51.79°N、123.08°E，海拔950 m。样芯带回实验室，经干燥、打磨等预处理后，选取22棵树的39根样芯用作重建。对选取的样芯，进行交叉定年，并获取宽度和密度数据，建立树轮晚材最大密度（MXD）年表。此外，我们从国家气象局气象数据共享网^[1]上收集了漠河站气象数据，并将树轮晚材最大密度年表与该站月平均气温和月降水量做相关分析，确定该点樟子松树轮晚材最大密度形成的限制因子。根据相关关系，建立年表与气候要素转换方程，重建该区过去的气温变化历史（图2）。

图2 数据集研发技术路线图

下载

4 数据结果与验证

4.1 数据集组成

基于树轮晚材最大密度重建的过去233年来大兴安岭北部7-8月平均气温数据集包括：（1）采样点的地理位置数据；（2）39根样芯的树轮晚材最大密度测量值统计特征（表3）；（3）树轮MXD标准年表（图3）；（4）重建的1781-2013年7-8月的平均气温，及11年滑动平均值数据（图4）。

图3 树轮晚材最大密度年表图

下载

图4 重建的1781-2013年7-8月平均气温变化图

下载

4.2 数据结果及验证

树轮晚材最大密度标准年表时间长度292年（1722-2013年）。选取1901-2000年作为所有序列的公共区间进行分析发现，所有序列间平均相关系数为0.408，同一棵树不同序列间的平均相关系数为0.56，不同树之间序列的平均相关系数为0.405，信噪比（SNR）22.09，样本量对总体代表性（EPS）为0.957。这些数据与该区域其他树轮密度年表类似[12]。该分析结果表明此年表可以用来进行气候分析。当1781年，具有7个样芯的时候，子样本信号强度（SSS）超过了0.85。因此，将1781年作为重建的起始年。

基于树轮晚材最大密度年表与气候要素关系，我们重建了大兴安岭北部地区1781-2013年7-8月平均气温变化历史（图4），其转换方程为T_mean7-8=3.46MXD+13.5，该模型对器测时期（1959-2013年）7-8月平均气温的方差解释量为31.1%。"留一法"交叉检验结果显示：符号检验原始序列通过了99%的检验，误差缩减值（RE=0.256）和乘积平均值检验t值均较高，重建和观测值的相关系数为0.52（n=55，p<0.01）；分段检验（校准期分别为1959-1988和1984-2013）结果显示：虽然仅有1959-1988校准期ST通过了99%检验，但误差缩减值（RE分别为0.279和0.319）和有效系数（CE分别为0.278和0.315）均大于零，并且验证期（1989-2013和1959-1983）重建和观测值相关系数分别为0.59和0.58。说明该重建稳定可靠，可以反应区域气温变化。

5 讨论和总结

本文利用采自大兴安岭北部呼中国家自然保护区的樟子松树轮样芯，获取了树轮密度数据，建立了树轮晚材最大密度年表，重建了该区过去233年7-8月平均气温变化历史。该数据集年代为1781-2013年，时间分辨率为年。本研究在空间上增加了气温重建分布点，为了解过去气候变化、寻求气候变化规律和预测未来气候变化提供了基础资料。

作者分工：李明启对数据集的开发做了总体设计；李明启和兰宇采集和处理了树轮密度数据；李明启和兰宇设计了模型和算法；李明启做了数据验证；李明启撰写了数据论文等。

利益冲突声明：本研究不存在研究者以及与公开研究成果有关的利益冲突。

参考文献

• [1] Fritts, H. C. Tree Rings and Climate [M]. London: Academic Press, 1976.
• [2] Liang, E. Y., Shao, X. M., Qin, N. S. Tree-ring based summer temperature reconstruction for the source region of the Yangtze River on the Tibetan Plateau [J]. Global and Planetary Change, 2008, 61(3/4): 313-320.
• [3] Shao, X. M., Huang, L., Liu, H. B., et al. Reconstruction of precipitation variation from tree rings in recent 1000 years in Delingha, Qinghai [J]. Science in China Series D: Earth Sciences, 2005, 48(7): 939-949.
• [4] Cook, E. R., Anchukaitis, K. J., Buckley, B. M., et al. Asian monsoon failure and megadrought during the last millennium [J]. Science, 2010, 328(5977): 486-489.
• [5] Yu, J., Shah, S., Zhou, G., et al. Tree-ring-recorded drought variability in the Northern Daxing'anling Mountains of Northeastern China [J]. Forests, 2018, 9(11): 674. DOI: 10.3390/f9110674.
• [6] Hughes, M. K., Schweingruber, F. H., Cartwright, D., et al. July-August temperature at edinburgh between 1721 and 1975 from tree-ring density and width data [J]. Nature, 1984, 308(5957): 341-344.
• [7] Jacoby, G. C., Ivanciu, I. S., Ulan, L. D. A 263-year record of summer temperature for Northern Quebec reconstructed from tree-ring data and evidence of a major climatic shift in the early 1800s [J]. Palaeogeography Palaeoclimatology Palaeoecology, 1988, 64(1/2): 69-78.
• [8] Briffa, K. R., Jones, P. D., Bartholin, T. S., et al. Fennoscandian summers from AD 500 temperature changes on short and long timescales [J]. Climate Dynamics, 1992, 7: 111-119.
• [9] Li, M., Huang, L., Yin, Z. Y., et al. Temperature reconstruction and volcanic eruption signal from tree-ring width and maximum latewood density over the past 304 years in the southeastern Tibetan Plateau [J]. International Journal of Biometeorology, 2017, 61(11): 2021-2032.
• [10] 白卫国, 王祝雄. 论我国东北林区森林可持续经营[J]. 林业资源管理, 2008(6): 7.
• [11] 唐亮, 吴东立, 苗微等. 东北地区食物安全可持续发展战略研究[J]. 中国工程科学, 2019, 21(5): 19-27.
• [12] Sun, Y., Wang, L. L., Chen, J., et al. Reconstructing mean maximum temperatures of May-August from tree-ring maximum density in North Da Hinggan Mountains, China [J]. Chinese Science Bulletin, 2012, 57(16): 2007-2014.
• [13] 李明启, 兰宇. 基于樟子松树轮晚材最大密度重建的大兴安岭北部7-8月平均气温数据集（1781-2013）[J/DB/OL]. 全球变化数据仓储电子杂志, 2022. https://doi.org/10.3974/geodb.2022.04.02.V1. https://cstr.escience.org.cn/CSTR:20146.11.2022.04.02.V1.
• [14] 全球变化科学研究数据出版系统. 全球变化科学研究数据共享政策[OL]. https://doi.org/10.3974/dp.policy.2014.05 (2017年更新).
• [15] Cook, E. R., Briffa, K. R., Meko, D. M., et al. The segment length curse in long tree-ring chronology development for paleoclimatic studies [J]. Holocene, 1995, 5(2): 229-237.
• [16] Cook, E. R. A time series analysis approach to tree ring standardization [D]. Tucson: The University of Arizona, 1985.